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Abstract. We investigate the partially asymmetric exclusion process (PASEP) 
with open boundaries when the reverse hopping rate of particles q = — 1, using a 
representation of the PASEP algebra related to the al-Salam Chihara polynomials. 
When q = — 1 the representation is two-dimensional, which allows for straightforward 
calculation of the normalization, current and density. We note that these quantities 
behave in an a priori reasonable manner in spite of the apparently unphysical value of 
q as the input, a, and output, j3, rates are varied over the physical range of to 1. 

As is well known, another two dimensional representation exists when < q < 1 
and abq = 1, where a = (1 — q)//3 — 1 and b = (1 — q)//3 — 1, and we compare the 
behaviour at q = — 1 with this. An extension to generalized boundary conditions 
where particles may enter and exit at both ends is briefly outlined. We also note 
that a different representation related to the q-harmonic oscillator does not admit a 
straightforward truncation when q = — 1 and discuss why this is the case from the 
perspective of a lattice path interpretation of the PASEP normalization. 
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1. Introduction 

Analytically continuing physical parameters such as temperature and field to unphysical, 
even complex values, is often done in the study of equilibrium statistical mechanical 
models, for instance in determining the Lee- Yang [TJ [2] and Fisher [3j |2j zeroes 
for partition functions whose scaling properties provide an alternative approach to 
characterizing phase transitions. It has been less commonly practised for non- 
equilibrium models, though Blythe and Evans have considered "normalization zeroes" 
for the asymmetric exclusion process (ASEP) by taking complex injection or removal 
rates at the boundary and recovered the known phase diagram for the non-equilibrium 
steady states [5j[6] from the scaling properties of these zeroes. More recently, Sasamoto 
and Williams [7] showed that choosing a negative injection rate at one boundary could 
relate the finite and semi- infinite partially asymmetric exclusion process (PASEP) over 
the region in which the finite PASEP was defined. In [7] the interesting general question 
of whether one can still give a probabilistic or physical meaning to the corresponding 
stationary distribution if one makes one or more parameters negative or complex in 
finite Markov chain models, such as the PASEP, was posed. 

The PASEP with a reverse hopping rate (defined below) of q = —1 appears to 
be an example of a model with a probabilistic interpretation since we find that not 
only does the normalization remain positive but quantities such as the current and 
particle density continue to behave in a reasonable manner. An additional motivation 
for choosing q = — 1 comes from consideration of finite dimensional representations 
of the PASEP algebra. Such n-dimensional representations have been obtained when 
< q < 1 along lines in the non-equilibrium steady state phase diagram of the PASEP 
given by ahq n ~ x = 1, where a = (1 — q)/a — 1 and b = (1 — q)/fl — 1 are convenient 
parametrizations of the injection and extraction rates a, (3. However, they also exist 
when q n — 1. These representations have not been considered before in the context 
of the PASEP precisely because the hopping rate takes an unphysical value but the 
study of the XXZ spin chain, which is closely related to the PASEP, at q n = 1 has 
been a topic of interest for many years (8 11. The simplest non-trivial example of such 



finite representations in the PASEP is found when q = — 1, giving the two dimensional 
representation whose properties we investigate here. 

We adopt a low-tech approach throughout, diagonalizing the 2x2 matrices of the 
PASEP algebra and explicitly calculating the current and density from these, which 
allows straightforward direct comparisons between the two-dimensional representations 
for abq = 1 and q — — 1. We briefly outline the existence of similar finite dimensional 
representations for q n = 1 (and q = — 1 in particular) when the boundary conditions are 
generalized. We also discuss the lattice path interpretation of the generating function for 
the PASEP normalization and note that a different representation of the PASEP algebra 
related to the q-harmonic oscillator does not admit finite dimensional representations, at 
least in an obvious manner, because it is no longer possible to implement the boundary 
conditions. 
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2. Model Definition and Basic Properties 

We start with the standard setup for the PASEP with continuous time Markov dynamics 
on a finite lattice of length N with open boundaries. Specifying the various injection, 




Figure 1. A typical particle configuration with allowed moves and rates in the PASEP 
model. Particles are injected on the left at a rate a when a space is available and leave 
at a rate f3 on the right. Their internal hopping rates are 1 to the right and q to the 
left when spaces are available. 



extraction and hopping rates for particles as shown in Fig. (1) defines the model. The 
model is clearly out of equilibrium since it supports a current driven by the injection 
and extraction of particles but non-equilibrium steady states can exist. The matrix 
product ansatz solution for the PASEP, 12 14 calculates the probability P(C) of a 
given configuration C in a non-equilibrium steady state by employing two vectors (W\ 
and \V) to sandwich a word X1X2 . . . X^ generated from an alphabet of two matrices 
D, E representing particles and holes in the configuration C. This gives 

= wx lX , x H \v) (21) 

Zn 

where the normalization for lattice of length N is 

Z N = (W\(D + E) N \V) = (W\C N \V) , (2.2) 

and we have defined C = D + E. These expressions are then related to the dynamics of 
the PASEP by imposing a quadratic algebra and boundary conditions on the matrices 
D, E and vectors (W|,|V) 

DE - qED = D + E , 

a(W\E = (W\ , (2.3) 
PD\V) = \V) • 

Exact solutions for any N may then be obtained by finding representations of the 
matrices and vectors or by normal-ordering the matrices to let them act on the 
appropriate boundary vector. The end result 13 is a phase diagram for the non- 
equilibrium steady states of the PASEP which contains a low-density (LD) a high 
density {HD) and a maximal current (MC) phase as the boundary rates a, f3 are varied, 
with the transition lines between the various phases being shown in Fig. (2). A novel 
feature is that a one-dimensional non-equilibrium system can exhibit (boundary driven) 
phase transitions, unlike its equilibrium counterparts in one dimension such as the Ising 
model. 
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Figure 2. The phase diagram of the PASEP. Here HD, LD and MC denote the high- 
density, low-density and maximal-current phases, respectively. The mean field line (a 
hyperbola) ab — 1 is indicated, below which no finite dimensional representations exist 
when < q < 1. We also show the dashed line abq = 1 along which the physical two 
dimensional representation may be defined, taking q = 1/2 for definiteness. 



In general, representations of the PASEP algebra are non-unique and will be infinite 
dimensional. The particular choice of representation which we find useful here is related 



to al-Salam Chihara [15] polynomials and is given by 16j-|18 
= (1,0,0,---) 



\V) 



1,0,0, 



and 



D 



1-q 



( 1 + 6 




V ; 

/ 1 + a 



E 



1-q 



1 + bq 




y/ci 1 + aq 





1 + bq 2 






\ 







y/02 1 + aq 2 



V 



/ 



where the various parameters which appear are defined as 

a = 1 , 

a 



c n = (l-q n )(l-abq 



(2.4) 



(2.5) 



(2.6) 
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With these choices the matrix C is 



C = D+E 







/ 2 + a + b y/c[ 

^fcx 2 + (a + b)q y/ci 

y/c~2 2 + (a + b)q 2 



(2.7) 



/ 



and the normalization Zjq in equ. (2.2) may then be calculated using various techniques. 
Similar expressions involving C, D and E, which we discuss below explicitly for both 
abq = 1 and q = — 1, may also be evaluated for the current and density. 

In the sequel we first reiterate the known results for the physical two-dimensional 
representation which exists when abq = 1, before going on to consider q — — 1. We 
compare the features of physical quantities such as the current and densities in both these 
representations and also discuss the representation related to the q-harmonic oscillator 
where, at least naively, continuation to q — — 1 does not appear to be possible. We 
conclude with some general comments on other representations with q n = 1 for n > 2 
and various possible generalizations. 



3. The physical two-dimensional representation: abq = 1, < q < 1 

The study of finite dimensional representations of the PASEP algebra and other 
quadratic algebras associated with different react ion- diffusion models pre-dated the 



exact matrix product solution of the PASEP 19 21 . It is useful to think of C as a 



transfer matrix (indeed, when the PASEP is mapped onto a lattice path model, as 



discussed in section 7, it is a transfer matrix) so looking at equ. (2.7) we can see that 
an n-dimensional representation will exist if c n = (1 — q")(l — abq 11 ^ 1 ) = since this 
decouples the matrix into blocks. Only the block in the upper left-hand corner will 
give a non-zero scalar product with the vectors (W\, \ V). One way this can occur is if 
abq™ -1 = 1, which restricts us to hyperbolic lines in the a, b plane of the PASEP phase 
diagram for physical values of q (i.e. < q < 1). The simplest of these is n = 1, giving 
ab = 1 which is just the mean field line for the PASEP, as indicated on Fig. (2). All the 
other finite dimensional representations may be defined on hyperbolae which lie above 
this, for instance the two-dimensional representation when abq = 1, which is shown in 
Fig. (2) for q = 1/2. 

It is instructive to compare the two-dimensional representation at q = — 1 with this 
other, physical, two-dimensional representation which exists along the line abq = 1 for 



< q < 1. This latter has been discussed by both Essler and Rittenberg 19 and Mallick 



and Sandow 20 in some generality, including the possibility of extraction at a rate 7 at 
the left boundary and injection at a rate 5 at the right boundary. The overall structure 
of the phase diagram is not altered by this embellishment (at least when q < 1) so we 
largely stick to the case of 7 = 5 = for simplicity. 

The line abq = 1 can cross the transition between the HD and LD phases in Fig. (2) 
as a, (3 are varied, so we must choose either a = 1/bq or b = 1/aq depending on the 
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dominant eigenvalue. Taking the former for definiteness (i.e. the HD phase), C is given 
by 

l+b 2 q+2bq J_ 

ri — ( bq(l-q) y/q \ (n 

° — I i i+fe 2 +2b I y 6 - 



and Z) and £7 by 



1+6 



n _ (1-9) 

^ 

V U (1-9) 

£ = ( l 7>* . (3.9) 

\ v 7 ? (l-<?) 6 / 

The vectors (W| and | V) are simply truncated to two components 

(W\ = (1,0) |V> = (l,0f. (3.10) 

For the LD phase we exchange b -H- a, and it is straightforward to verify that both 
the quadratic algebra relation and boundary conditions of equ. (2.3) are still satisfied in 
both phases. 

The off-diagonal elements of C in this representation are imaginary but the 
expressions contributing to physical quantities such as the current and densities are 
real and positive. The eigenvectors of C are 

Ai = (y + b + 2 



1 — q \b 

A 2 = t^— ( ^ + bq + 2 ) (3.11) 



1 — q \bq 

and it may thus be diagonalized using 



S= I VI ] (3.12) 



1 1 



and 



This gives 



= -i_ ib f r - 1 ] . (3.13) 
6 2 g - 1 \ -iby/q b q ' 



A = S~ l CS = I ^ ° I (3.14) 



A 



2 



and we can also define D and i£, which are both lower diagonal in the new basis 



19 



d = s-'ds I ■' v~i«j_ 



(1-9) 
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which, as we shall see below, simplifies the expressions for densities (and correlators) 
somewhat by comparison with q = — 1. The vectors (W\, \V) in the diagonalized basis 
are given by 



(W\ = (W\S= l-ijqb, 



where 



\V) =5- 1 |V)=m(l,-l) T 
ib^/q 



m 



(3.16) 



(3.17) 



b 2 q - 1 

and the combinations contributing to physical quantities again contrive to be real and 
positive. 

The first quantity of interest is the normalization, which is given by 

Z N = (W\C N \V) = (WISS'^SS-^V) = (W\A N \V) (3.18) 
and may be evaluated by using the definitions above to give 



Z N — 



b 2 q 



N 



b 2 q 



Ai 



N 



(3.19) 



In writing it in this form we have noted that Ai > A2, since b > 1/ ^/q along abq = 1 in 
the HD phase. We can thus define a correlation length 



= ln 



and write the normalization as 
Zn = ^1 



p — a exp I — 



N 



(3.20) 



(3.21) 



where 



b 2 q 



(Vq-l) 
1 



~ (6^-1) (3 - 22) 
which clearly disentangles the bulk asymptotic term from the finite length corrections. 
The current Jn is defined in the usual way as 

Zn-i 



Jn — 



Z 



N 



so asymptotically it is 



(3.23) 



(3.24) 



and we recover an expression for the current in the HD phase which agrees with the 
mean- field value. 
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The density may be calculated in a similar manner to give 
Z N ( Tj ) = (WIC^DC^IV) 
= (W\A?- 1 DA N - j \V) 
which may be written as 

Z N (rj ) = A^ -1 jx — to exp 



(N 



— C ex P 



(j - 1) 



(3.25) 



(3.26) 



where 



c 



b 2 q(l + b) 
(1-qWq-l) 

(1 + bq) 
(6»g-l)(l-g) 
6 



We can see from equ. ( |3.26[ ) that asymptotically in the density r starts at the value 

T = -r ^exp(-l/0 

pAi pAi 

at the left boundary and increases exponentially to 
\i b 
pAi ~ 1 + b 

in the bulk. We plot the density profile in the HD phase directly calculated from 



(3.28) 



(3.29) 



equ. (3.25) in Fig. (3), where we have taken b = 4, q = 1/2, iV = 100, which shows this 
behaviour clearly. We may evaluate two point correlators in a similar manner, giving 



0.785- 



10 20 30 40 50 60 70 80 90 

j 



Figure 3. The variation in the density (tj) along the lattice when abq — 1, for 
b — 4, q — 1/2, N — 100 (in the HD phase) showing the exponential increase to the 
bulk value of 4/5 from the left hand boundary. 
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exponential correction terms which continue to play a role only near the left boundary 
for the connected correlator. 



Similar expressions may be evaluated for the LD phase 19 , here the density stays 
low in the bulk before increasing exponentially at the right boundary and the connected 
correlator corrections are non-zero only close to the right hand boundary. The results 
might be summarized by saying that for both the HD and LD phases one of the 
boundary rates (b or a respectively) dominates the behaviour and corrections to the 
bulk are seen at one of the boundaries (left or right respectively) only. The maximal 
current phase is inaccessible to the two-dimensional representation along abq = 1 as can 
be seen in Fig. (2). 

4. Another two-dimensional representation: q = —1 

Since c n = (1 — q n )(l — abq 11 ^ 1 ) = when q n = 1, n-dimensional representations also 
exist if q is taken to be a root of unity. Let us ignore any qualms about the physical 
interpretation of such values of q and simply set q — — 1 to find another two-dimensional 
representation with the following expressions for C 

C " \\^2^M> l-*±* 1 (I " >U) 
and the individual D, E 

§(1 + 6) \y/2~- 



2V 1 ' 2 

|(1 




D = 

The expressions for the vectors remain the same as the other two-dimensional 
representation at abq = 1 



(W\ = (1,0) | V) = (1, 0) 



T 



and both the quadratic algebra relation and boundary conditions of equ. (2.3) can be 
seen to be satisfied here too. 

It is simpler to diagonalize C for explicit calculations just as it was for abq — 1. C 
now has the eigenvectors 

A+ = l + ly/2 + a? + b 2 

A_ = 1 - ^V2 + a 2 + b 2 (4.32) 
so we may use 

t t 



S = [ r -^- b a + b+r ) (4.33) 
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and 



2r 
_ J_ 

2r 



a+b—r 

2r 
i+b+r 
2r 



(4.34) 



to diagonalize, where we have defined r = y2 + a 2 + b 2 and t = y2 — 2ab for 
conciseness. The natural range of the parameters a, b when q = — 1 is 1 < a, b < oo, 
since the physical injection and extraction rates satisfy < a, < 1. The eigenvalues 
thus range over 2 < A + < oo, — oo < A_ < — 1. We would therefore not expect any 
phase transitions, which require A + = A_, since this is not possible for physical values 
of a, (3 unlike the abq = 1 representation. 
Diagonalizing, we have 





A = s^CS -- 
and similarly for D, E 

D = S^DS 

E = S^ES 



1 + 




1 - 



(4.35) 




1+6 2 




b+ l + b* 


r 

1+b 2 




r 

I 1+b 2 
r 


r 




1+a 2 
r 


b 


- r \ 


1+b 2 
r 


1 





(4.36) 



The vectors (W\, \ V) in this basis are 

(W\ = (W\S = \{a + b + r,a- h 



t 



t 



\v) =s- 1 \v) = -n,-i) 



2r 



(4.37) 



We can see that D and E in equ. (4.36) are no longer lower diagonal, unlike the D and 
E in equ. (3.15) when abq = 1, which will have consequences for the behaviour of the 
density and correlators. 

The first quantity of interest is again the normalization Z^f, which may be evaluated 
to give 

1 



2r 



(A- 



(a + b + r) - (a + b 



A_ 



N' 



We can define the correlation length in this case by 



1 

I 



In 



A, 



IA_| 



and write 



with 



A 



N 



-l^crexp 



N 

7 



(4.38) 



(4.39) 



(4.40) 



a 



a + b + r _ a + b+V2 + a 2 + b 2 

2~r " 2^2 + a 2 + b 2= 
a + b-r a + b- V2 + a 2 + b 2 
~ 2V2 + a? + b 2 



2r 



(4.41) 
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Although A_ is negative, Zn itself is positive for all N (both even and odd) as can be 
seen by expanding equ.(4.38). For instance, Z 3 is given by 

Z 3 = 5/2 + (7/4)a + (7/4)6 + (3/4)a 2 + (3/4)6 2 



+ (l/8)a 3 + (l/8)a 2 6 + (l/8)a& 2 + (l/8)6 3 . 

This in turn means that the current 

Zn-i 



(4.42) 



J 



N 



is a positive quantity. As we can see in Fig. (4) it is a smoothly decreasing function of 




Figure 4. The current Jn for N = 100, which is typical, plotted for a, b = 1 . . . 10 



a, b dropping from a value ofJ=l/2ata = 6= l (i.e. a = (3 = 1) to zero at large a, b. 
Asymptotically in iV it is given by 

J = Jn^oo = — = , 1 . (4.43) 

Remembering that a = 2/a — l, b = 2/(3 — 1 when q — — 1, increasing a, b from a = b = 1 
corresponds to decreasing the injection and extraction rates a, (3 so the observed falloff 
of J and the values it takes in Fig. (4) are not physically unreasonable. 

When q = — 1 the current J in equ. ( |4.43 ) depends on both a and b, which is not the 
case when < q < 1. For the latter it is given by (1 — q)a/(l + a) 2 and (1 — q)bj (1 + b) 2 
in the low and high density phases respectively and by J = (1 — g)/4 in the maximal 
current phase. A naive continuation to q = — 1 from the maximal current phase which 
extends to a = {3 = 1 gives the value of J = 1/2 observed here for q = — 1 at a = {3 = 1. 
Zn and hence Jn contain exponential corrections which oscillate in sign, a feature which 
is also apparent in the density, which we turn to next. 

If we evaluate the density 

. . (W\Ci- l DC N -i\V) 
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(wis-^sc^s-^iSDS-^isc^^s-^siv) 



J N 



(W\A j - 1 DA N ^\V) 



(4.44) 



the individual terms in the numerator no longer remain positive, for example 

(W\CDC\V) = 5/4 + b + (3/4)o + (5/8)6 2 

+ (l/8)a 2 -(l/8)a 2 b + (l/4)a6 2 + (l/8)b 3 
but the numerator as a whole does stay positive. In general 
Z N { Tj ) = {WIC^DC^IV) 



(4.45) 



(a + 6 + r)A 
4r 



AT-l 



1 + 



(a + 6 - r)A 
4r 



iV-l 



1 + 6 2 



1 + 6 2 



(4.46) 



which, if we adopt a similar notation to equ. (4.40), may be written as 

(N-jy 



Z N ( Tj ) = X 



N-l 



l) N - j 9exp 



(4.47) 



l^exp 



(iV-l) 



Cexp 



(j ~ 1) 
I 



This should be contrasted with the corresponding expression for abq = 1 in equ. (3.26) 



which contains one less exponential term because of the lower diagonal form of D in 
that case. The various coefficients in the above are given by 

(a + b + V2 + a 2 + b 2 ) (_ 1 + b 2 



4^2 + a 2 + b 2 



1 



c 



{a + b + V2 + a 2 + b 2 ) 

4^2 + a 2 + 6 2= 
(a + 6 - V2 + a 2 + b 2 ) 



(a + b 



4^2 + a 2 + b 2 
V2 + a 2 + 



b 2 ) 



4^2 + a 2 + b 2 



b + 



V2 + a 2 + b 2 

l + b 2 
V2 + a 2 + 6 2 

1 + & 2 
V2 + a 2 + 6 2 

1 + & 2 
V2 + a 2 + 6 2 



(4.48) 



As with the current, the behaviour of the central (bulk) density appears to be physically 



reasonable. From equ. (4.47) and Fig. (5) we can see that the central density (rso) 



depends on both a and 6, rather than just b as in the HD phase and a in the LD phase 
when < q < 1. From Fig. (5) it is also apparent that (T50) decreases monotonically with 
increasing a (i.e. decreasing injection rate) and increases monotonically with increasing 
b (decreasing extraction rate). The bulk value of the density is given asymptotically by 

T = J- = l + b 2 + V2 + a 2 + b 2 

pA+ V2 + a 2 + b 2 (2 + V2 + a 2 + b 2 ) 
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Figure 5. The central value of (tj), for j = 50, N = 101, calculated directly from 



cqu. (4.44) and plotted for a, b = 1 



10. 



so, in particular, we can see that r = 1/2 when a = b, as evidenced by the central line 
along the surface in Fig. (5) and in the interior values in Fig (6). The asymptotic formula 



of equ.(4.49) is already an extremely good match to the directly calculated values shown 



in Fig. (5) when N = 101 over the full range of a, b. Note that, as is clear from Fig (5) 



and equ. (4.49), r is not symmetric in a, b, unlike the current, but the density of holes 
1 — r is given by exchanging a <H- b 

V2 + a 2 + b 2 (2 + V2 + a 2 + b 2 ) 
If we now look at the density profile along the lattice for given a, b there are 
oscillating positive and negative exponential corrections at both ends, when j ~ 1 and 
j ~ N, in contrast to the two-dimensional representation for abq = 1 



19,20 where 



corrections to the bulk are apparent at only one boundary. The oscillating signs when 
q = —1 in the correction terms to the bulk density will still be present even if we consider 
only odd or even length lattices. A similar story holds for two point correlators, which 
continue to display "two-sided" oscillating sign corrections, in contrast to the corrections 
in the HD and LD phases when abq = 1 discussed in the previous section. 



5. Finite dimensional representations when 7,5 ^ 

The PASEP boundary conditions may be generalized to allow wrong direction injection 
and removal rates. Particles are now injected at a rate a and removed at a rate 7 at the 
left boundary and particles are removed at a rate (3 and injected at a rate 5 at the right 



boundary. Tridiagonal representations of D and E with non-zero 7 and 5 still exist 17 
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10 20 30 40 50 60 70 80 90 



Figure 6. The variation in the density (tj) along the lattice when q = — 1, for 
a = 20, b = 20, N = 100, showing the oscillating sign exponential corrections at both 
ends. The large a, b values were chosen to make the oscillating corrections clearly 
visible. 



and they may be written as 

i 



D 



1-q 



4 







1 + d\ d\ 
d\ 1 + d\ 



E 



1-q 



1 1 + 4 



V 



l + e\ 



el 



1 + e 



/ 



(W\ = (1,0,0, ■■■), 

where d\, d\, 

e n> e n are given by 



\V) = (1,0,0, ■■■)' 



(5.51) 



(5.52) 



di 



„n-l 



(1 _ q *n-*abcd){\ _ q 2n abcd j 

x [6d(a + c) + (6 + d)g - a&ceZ(6 + d)q n ~ 1 - {bd(a + c) + abcd(b + d)}q n 
- bd(a + c)q n+l + a& 2 cd 2 (a + c)^™" 1 + a6cd(6 + d)g 2n ] , 



,n-l 



(1 - g 2 "- 2 a6cd)(l - q 2n abcd) 

x [ac(6 + d) + (a + c)<? - abcd(a + c)g™ _1 - {ac(b + d) + abcd(a + c)}g r 
- ac(5 + ci)g n+1 + a 2 6c 2 d(6 + d)^ 1 + abcd{a + c)g 2n ] , 



di = 



1 

1 — q n ac 



A 



e" = 



q n ac 
1 — q n ac 



A 



(5.53) 
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q n bd 

1 _ q n M An ' 
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1 A 

1 - q^bd n ' 



which involve the further parameters 

1 



d 
A 



2a 
1 

1 

2a 
1 

2~j3 



(1 - q — a + 7) + _ 1 ~ a + 7) 2 + 4cry 
(1 - g - f3 + 5) + v / (l-g-/3 + 5) 2 + 4/35 

(1 - g - a + 7) - a/(1 - g - a + 7) 2 + 4c*7 
(1 - g - /3 + 5) - v/(l-g-/3 + 5) 2 + 4/35 



(5.54) 



(l-g n - 1 a6cd)(l-g"+ 1 )(l- g "afc)(l-g n ac)(l-g' 1 ad)(l-g"fec)(l- g "fed)(l-q"cd) 
(l-g 2ll - 1 afecd)(l-g 2 "afecd) 2 (l- ( j 2 ' l + 1 afecd) 



■,1/2 



In this case the polynomials associated with the representation are the Askey- Wilson 
polynomials [22], of which the al-Salam Chihara polynomials are a specialization. 

We note that when 7 = 5 = 0, a and b revert to the definitions of previous sections 
and c = d = 0. Although the expressions above are considerably more complicated 
than the 7 = 5 = case the structure of C, D and E is still similar. The diagonal and 
off-diagonal elements of C are now given by d n and a/£^ where 
2 + di+ei 



1-? ' 

dl + el){dj + ej) 
(1-g) 2 



(5.55) 



Explicitly, 



C n +1 



1 - gn^otedJCl - q n+1 ){l - q n ab){l - q n ac)(l - q n ad){\ - q n bc)(l - q n bd)(l - q n cd) 



(1 - q 2n - l abcd){\ - q 2n abcd) 2 {\ - q 2n+l abcd){\ 



(5.56) 



and n-dimensional representations will still exist if one of the terms in the numerator 
of c n is zero. This may occur in various ways: abq n ~ x = 1 as before; variations thereof 
involving at least one of c, d such as acq n ~ l = 1 or cdq n ~ x = 1; or even a four-parameter 
condition abcdq n ~ 2 = 1. The latter two possibilities are not available when 7 = 5 = 0. 
We can also see that n-dimensional representations will continue to exist when q n = 1, 
since a factor of 1 — q n remains present in the numerator of c n . 

Curiously, conditions of the form abcdq n = 1 and q n = 1 are both encountered when 



solving the open PASEP using the Bethe ansatz 23 27 , but the exact relation between 



finite dimensional representations of the matrix ansatz (indeed, the matrix ansatz in 
general) and the Bethe ansatz solutions is unclear. 



6. A representation that does not truncate gracefully at q— — 1. 



An alternative representation the PASEP algebra related to the q-harmonic oscillator 
28 and q-Hermite polynomials was employed in the original solution of the PASEP 
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in 



13 . In this D and E are written as 
1 



D 



E 



, +1 ( 

1 — q 1 — g 

1 1 . 



t 



1 — q I — q 

where a, a) satisfy q-boson commutation relations as a consequence of 
algebra 

— ga^a = 1 — q 



aa' 



a' \n 



)= (l-gn+^ln+l) 



a| n ) = (l-g n ) 1/2 |n- 1) 
a|0) = . 

The matrices D and are given explicitly in this representation by 



D 



l-q 



1 








\ 



'l-q 2 
1 



E 



1 



V 



l-q 

and the vectors {W\, \ V) by 
(W'ln) = re 



V : 

/ 1 






1 








(n\V) = re- 



V(T,q)n 



with 



V (q; q)n 

1 _ (ab) n 



n=0 



We have used the standard notation for (shifted) g-factorials in the above 



n-l 



(a;q)n = n^ 1 ~ aqJ 

3=0 

The boundary conditions 

a(W\E = (W\ 
0D\V) = \V) 
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are implemented in a slightly different manner in this representation compared with the 
previous section. We can see that the action of D on \V), for instance, 



D\V) 



1-g 



/ 1 







1 








T -q 1 



\ 



J 



( 



\ 



fY^q 
b 2 



\/(l-9)(l-? 2 ) 



\ 



(6.64) 



is to shift up a component of \V) and add it to the component above in order to satisfy 
PD\V) = \V). 

Although the g-boson relations collapse to suitably scaled fermionic anti- 
commutation relations for a two-state system when q — — 1 



dd* 


+ d f d = 2 


0*10) 


= V2|1> 


d|i) 


= V2\0) 


d|0) 


= 


d f |i) 


= 



(6.65) 

and the matrices C,D and E decouple into 2x2 blocks, the vectors, for instance \V) 
in equ. (6.64), and normalization n can be seen to be ill-defined when q = —1. This 
particular representation is therefore not suitable for continuing to q — — 1 (nor, indeed, 
q = 1 [l4]). In the next section, where the relation between the PASEP and lattice path 



models is outlined, we see heuristically why this is so. 



7. The PASEP and Lattice Paths 

The generating function of 

Z{z) = Y,Z N z N 

N 



(7.66) 



can be thought of as a "grand-canonical" normalization where the PASEP 
normalizations for various iV are combined with a fugacity z. The closest singularities 
of Z(z) to the origin can then be used to determine the asymptotic behaviour of Z^, a 



common technique in analytic combinatorics 29 , 30 



It is also useful to interpret 2(z) as the generating function for a model of weighted 
which gives some insight into the PASEP phase diagram. From 



lattice paths 31 36 



this point of view the tridiagonal matrix C in the particular representation of equ. (2.7) 
is a transfer matrix for Motzkin paths composed of diagonal up and down steps with 
with weights yfc^zjiX — q) at height n and horizontal steps with weight d n z/{l — q) 
at height n. A path of length n will be weighted by z n and some path-dependent 
product of c n and d n , so by tuning these via a, (3 and q we can change the dominant 
paths contributing to the ensemble. Since only the first components of (W\ and \V) 
are non-zero the paths start and end on the horizontal axis such as the path shown in 
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Figure 7. A Motzkin path of length 12 which reaches a height 3. The horizontal axis 
is shown dashed so the horizontal step on the axis is visible and the initial and final 
point arc indicated 



Fig. (7). The horizontal steps may be composed into two different "colours", weighted 
by (1 + aq n )z/{\ — q) and (1 + bq n )z/(l — q) respectively, so the model is really one 
of bi-coloured Motzkin paths, which are of considerable interest to combinatorialists 
because may be put into bijection with numerous other objects. 

In the lattice path model the dominant paths contributing to the normalization 
in the HD and LD phases contain mostly one colour of step bound closely to the 
horizontal axis, where the binding energies are given in terms of a and 0. The maximal 
current phase, on the other hand, corresponds to unbound paths in which entropy 
dominates 14 . Restricting C to be finite in this context means that the heights of the 
paths in the ensemble are restricted to lie below a ceiling determined by the dimension 
of the representation. It is therefore no surprise that finite dimensional representations, 
which cannot describe unbound paths, do not see the maximal current phase. 

The generating function of in the q-harmonic oscillator representation could also 
be viewed as a Motzkin lattice path model since the matrices there are still tridiagonal, 
but in this case the form of the vectors in equ. (6.60) means that the paths may start 
and end at any height, showing why a consistent truncation to height zero and one paths 
does not appear to be possible in that case. 



8. Conclusions 



We have seen that a two-dimensional representation of the PASEP quadratic algebra 
at q = — 1 presents broadly reasonable physical behaviour, although the oscillating sign 
exponential corrections at both ends are unusual. We compared its properties in detail 
with the two-dimensional representation which exists when abq = 1 and < q < 1. 
The behaviour when the boundary conditions were extended to non-zero 7 and 5 was 
found to be similar, in that two-dimensional physical representations existed when 
parameters were restricted appropriately, e.g. abq = 1 or abed = 1, alongside another 
two-dimensional representation when q — — 1. 

We found that a naive continuation to q = — 1 was not possible for the 
representation of the PASEP related to the q-harmonic oscillator because of the form 
of the vectors (W|,|V) in that case and presented some heuristic arguments from a 
weighted lattice path interpretation of the PASEP as to why such behaviour might have 
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been expected. 

Exploring n-dimensional representations for q n = 1, n > 2 would be a possible 
extension of the current work, but we note that in such cases the Zn is no longer 
necessarily real when expressed as a function of real a, {3 so extracting a physical 
interpretation for currents, densities and correlators may not be quite so straightforward. 
Having already allowed the unphysical value of q = —1 in the two-dimensional 
representation, one might also take unphysical values of a, b in that case. It would then 
be possible to arrange A+ = A_ = 1, and hence a phase transition, at a = b = i. Another 
avenue which might be pursued is the use of larger representations with g n = 1, n > 2 
to approach the symmetric limit of the PASEP at q = 1 through large finite dimensional 
representations of the PASEP algebra. 

A further task is, of course, finding a plausible physical setting for the two- 
dimensional q = —1 representation of the PASEP discussed here. We have been unable 
to do so ourselves thus far. 
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